Least action principle for envelope functions in abrupt heterostructures 
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We apply the envelope function approach to abrupt heterostructures starting with the least ac- 
tion principle for the microscopic wave function. The interface is treated nonperturbatively, and 
our approach is applicable to mismatched heterostructure. We obtain the interface connection rules 
for the multiband envelope function and the short-range interface terms which consist of two phys- 
ically distinct contributions. The first one depends only on the structure of the interface, and the 
second one is completely determined by the bulk parameters. We discover new structure inversion 
asymmetry terms and new magnetic energy terms important in spintronic applications. 
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The envelope function method is a powerful tool which 
has been widely used to describe and predict various 
effects in semiconductors. It is normally applicable to 
materials with translation invariance (allowing for the 
expansion of the wave function into Bloch functions) 
and to slowly varying potentials. There are two com- 
peting approaches to extending this method to abrupt 
heterostructures 0] taking into account interface-related 
effects. The first one is to impose appropriate boundary 
conditions (interface connection rules) on the envelope 
wave function at the interface 0. El [J Q . Another possi- 
bility is deriving the exact envelope function differential 
equations which are valid near the interface and which 
contain the iterface-related terms [f| ■ The second ap- 
proach is more detailed, and it requires a lot more infor- 
mation on the microscopic structure of the interface. Up 
to now, it has only been applied to the lattice-matched 
heterostructures where the interface is a weak perturba- 
tion. In this case, it has been shown that connection 
rules and differential equations are equally valid repre- 
sentations of the interface behavior. 

It is the aim of this letter to present an extension of the 
envelope function method which treats the interface non- 
perturbatively, and which is applicable to mismatched 
heterostructures. It turns out that the best approach to 
the problem is via the Lagrangian variational principle 
which encodes the Schrodinger equation. The advan- 
tage of this method is that both the Hamiltonian and 
the boundary conditions at the interface are contained in 
the averaged variational functional. The resulting k ■ p 
heterostructure Hamiltonian coincides with the ordinary 
k ■ p Hamiltonians on two sides of the interface. In addi- 
tion, it contains short-range interface (SRI) terms which 
are the main object of our study. We show that the 
SRI terms consist of two physically distinct contribu- 
tions. The first one is represented by the Hermitian in- 
terface matrix. Its components are directly connected to 
the parameters of boundary conditions for the envelope 
functions 0] and determined by the microscopic structure 



of the interface. The second contribution is completely 
determined by the bulk parameters of the materials. It 
includes new structure inversion asymmetry (SIA) terms 
and new SRI magnetic terms that are additional to the 
well known Rashba SIA terms and to the macroscopic 
magnetic terms, respectively. Taking them into account 
is important for various mechanisms of spin polarization, 
spin filtering and spin manipulation that play a key role 
in semiconductor spintronic applications 

In this letter we consider a model of a semiconductor 
heterostructure made of two homogeneous semiconductor 
layers A and B of characteristic length L. The layers are 
joined by a thin boundary region II of the width d w 
flo <C L (see Fig. [IJ, where ao is the lattice constant. 
We work in the single electron approximation, and we 
denote by U (r) the effective potential for electrons. U (r) 
coincides with periodic crystal potentials Ua,b{t) inside 
the bulk-like regions A and B, respectively. 

We start with the microscopic Lagrangian variational 
principle which encodes the stationary single electron 
Schrodinger equation. The corresponding Lagrangian 
density is of the form, 



£($*,$) = (E- J7(r))|$(r)| 5 
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Here mo is the free electron mass and $ is the micro- 
scopic spinor wave function. To simplify the presenta- 
tion we first neglect the spin-orbit terms in Eq. Q . The 
variational principle reads, 



SS = S [ d a r£($*,$) = 0. 



(2) 



where variations 6$ and 8&* are independent of each 
other and vanish at the outer boundaries of the in- 
tegration region SI = A + H + B. The variational 
principle implies the microscopic Schrodinger equation 
Hmicro&ir) = E$(r) with the microscopic Hamiltonian 
Hraicro = (p 2 /2mo + U(r)), where p = —ih'V is the mo- 
mentum operator. The microscopic probability flux den- 
sity j = ($*(r)p<I>(r) +p*$*(r)<I>(r))/2mo is conserved: 
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ties Ca,b are obtained from the microscopic Lagrangian 
by averaging over the unit cells in A and 7?, respectively. 
They have the form: 



£(**,*) = £|1<| 2 - **C* - — V^A^V,,* (3) 
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Here A^" = D>* v +£)"», K*> = e mv D» v , r],fj,,u = x,y,z 



and 



is Levi-Civita anti-symmetric tensor. The sur- 



face Lagrangian is nonlocal and it can be written as 
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FIG. 1: Sketch of the planar heterointerface between A and 
B semiconductor layers. II denotes the boundary region. 
Ua(Ub) is the crystal potential in A(B). 



Vj = 0, the microscopic wave function 4> is continuous 
everywhere in the heterostructure. 

It is our aim to pass from the description in terms of 
the microscopic wave function $ to the envelope func- 
tion approximation. To this end, we use expansions 
^ A ' B (r)u A ' B within A and B. Here 



*(r) = Zn=f 

u A ' B , n = Na.b are the Bloch functions at extremum 
points of the bulk energy band structure. The Na,b 
component envelope functions ty A ' B (r) are smooth in the 
A and B regions where they satisfy matrix Schrodinger 
equations H A > B (k)V A ' B (r) = E^ A > B (r). Here H A > B = 

are the standard fc • p Hamilto- 



C- 

nians including the terms up to the second order in the 
wave vector operator k = — iV. The matrices C, 73^ and 
v (/i ,v = x,y, z) are Hermitian Na.b x Na,b tensors 
of rank 0,1, and 2, respectively. The Hamiltonians H A ' B 
give a direct description of the Na,b basic bands as well 
as the contributions of the remote bands in the second 
order of perturbation theory 9]. Note that symmetry 
of the materials, the number of basic bands in the k ■ p 
approximation can be different on two sides of the inter- 
face. Moreover, parameters of bulk Hamiltonians H A ' B 
can vary significantly across the interface, so as it cannot 
be treated as a weak perturbation of the bulk problem. 

We fix once and for all the basic functions u A ' B in A 
and B, and we derive the h p version of Lagrangian vari- 
ational principle. The result has the form SS^ P = SSa + 
8Sb + SS SU1 - — 0, and it contains N = Na + Nb indepen- 
dent variations of the envelope wave functions S^ A , S^ 13 . 
Here SS a ,b = f A>B d 3 r£(<P*, $) = j AB d 3 r£ AtB {**, *) 
and SS sur — S J n d 3 rS(x)£ sur , where S(x) is the Dirac 
delta-function. The bulk multiband Lagrangian densi- 



with p = (y, z), and d = a + b 



where 

~ _ / y A (p,x~a) 

{* B (p, X + b) 
(see Fig. [JJ. The energy independent hermitian N x N 
interface matrix T sur depends on the symmetry of both 
bulk materials and of the interface. It can be constructed 
by using the method of the invariants or calculated di- 
rectly via the microscopic modeling of the potential U(r) 
in the interface region n (the details will be presented 
elsewhere) . 

The effective Lagrangians La.b together with £ sur con- 
tain all the relevant information about the bulk and 
interface properties of the heterostructure. Applica- 
tion of the least action principle SS^p = generates 
the Schrodinger equation Hab^ — with the com- 
plete heterostructure k ■ p Hamiltonian Hab and the 
general boundary conditions (GBC) to be imposed on 
<F a = * A (p, -a) and W fc = * s (p,6). The GBC (see 

iV T ^ a \ = 2T^ / * Q 

-iV T v b J h V * b 

iV T ^ a ) = TtI { iV T ^ b ) ' Where thC com P° nents 
of the 2~N a x 2Nb transfer matrix T tr (see 013) can be 
readily expressed via the components of the surface ma- 
trix T sur (see also 01)- Here r is the normal vector to the 
interface, V T = r -V, and the envelope velocity operator 
Vnm^m = 2i/%(dC/d V** ) can be written explicitly as 



[2j) can be written as 



V = B 



dk 
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The last term is new in comparison to 0]. The corre- 
sponding extra term in the envelope flux density J{r) — 

1/2 (**K m * m + V n (V nm V m )*) is proportional to V x 



(ty*K^) and does not alter the continuity property 
V • J = 0. It is straightforward to verify that J" 13 — 
t ■ J Q/3 (r) = const, where a and (3 label two functions 
^ a and ^ satisfying the same GBC (see Ref. 2]). This 
ensures that the heterostructure k ■ p Hamiltonian 77^ 
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is self-adjoint. It has the form: 
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where H A = Q{-x)H A (p, x - a), H,f = 0(x)H A {p, : 
b), Q(x) is the Heaviside step function and 
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We see that there are two physically distinct contribu- 
tions to the short-range interface (SRI) terms of the 
Hamiltonian Hab ■ The first one arises from the nonlocal 
surface Lagrangian £ sur and it depends on the properties 
of the interface via the energy independent parameters of 
the GBC. The other contribution comes from the veloc- 
ity operator V T . It is entirely determined by the bulk 
parameters and it arises from the nonvanishing variation 
of the bulk Lagrangians Ca,b at the interface. The im- 
portant feature of this contribution is the presence of 
the asymmetric K term. In homogeneous semiconduc- 
tors the asymmetric K term does not contribute to H 
in the absence of external fields (see ^3). Examples be- 
low demonstrate that the K terms in the Lagrangian of 
Eq. (PJ) and in the velocity operator of Eq. J^J become 
important if the symmetry is broken due to the presence 
of external fields or asymmetric interfaces. To emphasize 
this point we neglect effects of bulk inversion asymmetry. 

As a first example, we consider the effective mass 
Hamiltonian H(T 6 ) = E c + h 2 k 2 /2m for the spmor en- 
velope function ^ a (a = ±1/2), where E c is the bottom 
of the conduction band and m is the effective mass. Fol- 
lowing our method we introduce the effective mass La- 
grangian density 



£= {E-E c )\^\ 2 - — |V*(r)| 2 + £ S iA, 
2m 

which contains the asymmetric term 



£siA(r6) 
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obtained with K(Tq) = —(ig/2mo)cr, where a x , a y , a z 
are Pauli matrices, and g = go — g is the difference 
between free electron and effective electron g factors. 
Note that j ^ only if the spin-orbit splitting A of 
the valence band or A c of the remote conduction band 
is taken into account. We discover now that it is this 
asymmetric term £siA(r6) (missing in Refs. 0, 0|) 
that induces the spin dependence of the velocity operator 
V = (h/m)k — (ift,g/2mo)[crxk] and thus the spin depen- 
dence of the standard boundary conditions W = const, 



V T ^f = const at the interface (see 0,^3)- The short- 
range interface SIA term in the heterostructure Hamilto- 
nian Hab of Eq. © also results from ^sia^)- More- 
over, exactly this term £giA(r 6 ) produces the magnetic 
energy term —l/2/j,Bg{(TH) additional to the free elec- 
tron magnetic energy l/2/isgo( cr H) in the bulk semi- 
conductor in external magnetic field H. Here /is is the 
Bohr magneton. Next, it can be shown that the macro- 
scopic SIA term H so = aR[erxk)T, postulated by Rashba 
jlfij for the asymmetric 2D structure, is generated by the 
term £sia(F6). For this the dependence of g on the po- 
tential V = — \e\£x should be taken into consideration, 
where the average electric field E = £t characterizes the 
macroscopic asymmetry. In the eight band model for cu- 
bic semiconductors g = g — g r — 2E p A/3(E g — V)(E g — 
V + A), where E p is Kane energy, E g is a band gap and 
g r is a correction from remote bands, and the effective 
Rashba constant is an oc dg/dx\ x _ oc A. Using the 
expression for g in the 14 band model one finds that the 
correction to cxr is proportional to A c . 

Another useful example is provided by the degenerate 
valence band at the T point described by the envelope 
Hamiltonians obtained in |ll|. We consider two cases 
of A = and A — > oo. The remarkable property of 
the respective envelope Lagrangians obtained according 
to Eq. J3J is the existence of the asymmetric term with 
K(Ti5) = — i(l + 3k)/toqJ even in the case A = 0: 



-SIA 
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(l + 3re)V**[I x V* c 
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Here K is the magnetic Luttinger constant, I is the in- 
ternal angular momentum operator (/ = 1) and ^S a , 
a = 0, ±1, is the 3 component envelope function. The 
SIA component of the velocity operator V so = ifr(l + 
3n)/mo[I x fe] induces a new short range SIA term in 
the heterostructure Hamiltonian for the Ti5 holes as well 
as the /-dependent boundary conditions. This leads to 
the splitting of the heavy hole subband in asymmetric 
structures mediated by the interaction with light hole 
states. Note that it is this asymmetric term £sia(Fi5) 
which is responsible for the magnetic energy term oc 
/Xb(1 + 3k)(IH) in the bulk Hamiltonian of Ref. [ll|. 

In the case of A — * oo, the top of valence band is 
four- fold degenerate corresponding to the J = 3/2 sub- 
space of the total internal momentum J = J+l/2<r. We 
obtain the asymmetric term in the envelope Lagrangian 
with K(T$) = —i(2/3 + 2n)/m J — iq/moF, where q 
is cubically anisotropic magnetic Luttinger constant and 
F=(F a ,F v ,F t ) = (J*,J*,jS): 



-vw* [K(r 8 ) x v* a ] 



(11) 



Here ^> a , a = ±3/2, ±1/2, is the 4 component envelope 
wave function. The SIA component of the velocity op- 
erator V so = ih/m ((2/3 + 2k)[J x fe] + q[F x fe]) pro- 
duces a new short range SIA term in the heterostructure 
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Hamiltonian as well as the asymmetric contribution to 
the boundary conditions of the T 8 holes (see 0, llijh 
The very same asymmetric term £giA(r 8 ) induces the 
magnetic energy terms oc (JH) and oc q(FH) in the 
bulk Hamiltonian of Ref. as well as the macro- 

scopic SIA term H$° = BjJJ x k]E + f3 2 [F x k]E postu- 
lated recently in Ref. 17]. The cubically anisotropic 
constant q and, consequently, /?2 are proportional to 
A c and usually small. Considering the dependence of 
K = K r + E p /6(Eg + V) on V = —\e\£x, where K>t IS the 
contribution from remote bands, we are able to derive 
for the first time the effective SIA constant /3i for the r§ 
valence band: 
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We found that in contrast to the Rashba constant, /3i is 
not proportional to the spin-orbit splittings A or A c . 

The short-range SIA contribution to the zero field 
splitting becomes more pronounced in the case of the 
interface between very dissimilar materials. The general 
Hamiltonian Hab of Eqs. (|6I8|) enables to describe, for 
example, the interface coupling between Tq electrons and 
r 8 holes and the SIA effects in type II and type III quan- 
tum wells, and to take into consideration the microscopic 
asymmetry caused by the nonequivalence of two opposite 
interfaces (see |l8)h 

Finally, the developed approach reveals a new short- 
range magnetic energy term in the heterostructure 
Hamiltonian Hab m the presence of the in-plane exter- 
nal magnetic field. Indeed, choosing H || z we obtain 
the new term in the velocity operators V T ^ proportional 
to Hx(K z ^). In small fields this term can be treated 
as a perturbation. The short-range contribution to the 
Zeeman energy is proportional to the sum of the disconti- 
nuities of (^>*K Z ^>) at the interfaces. Unlike the zero field 
splitting, this interface magnetic energy contribution is 
present even in completely symmetric 2D structures as 
well as in spherical dots (see [lflj]V 

In conclusion, the variational least action principle al- 
lows to consistently extend the envelope function ap- 
proach to the heterostructures with abrupt interfaces. 
The SRI terms in the heterostructure Hamiltonian Hab 
and the GBC are equally valid representation of the in- 
terface properties and can be written for any interface be- 
tween dissimilar materials including the case Na ^ Nb- 
For lattice-matched heterostructures (Na — Nb = N) 
the obtained Hab allows direct comparison with previ- 
ously derived Hamiltonians 0, Q, 0] . The discretization 
of Hab for numerical calculations is straightforward and 
requires no additional symmetrization. All macroscopic 
and short-range interface SIA terms as well as the mag- 
netic energy terms in Hab originate from the asymmetric 
term £sia m the bulk envelope Lagrangian. 
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